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Abstract 

We compute various processes involving neutrinos in the initial 
and/or final state and we assume that neutrinos have energy momen- 
tum relation with a general power law E'^ = -\- ^^p" correction due 
to Lorentz invariance violation. We find that for n > 2 the bounds 
on from direct time of flight measurement are much more stringent 
than from constraining the neutrino Cerenkov decay process. 



1 Introduction 



The OPERA observation [1] of superluminal neutrinos has been ruled out by 
the measurement by ICARUS [2] which puts an upper bound on the super- 
luminahty of neutrinos at 5 < 2.3 x 10~^. An phenomenological consequence 
of neutrinos being superluminal was the observation by Cohen and Glashow 
(GC)[3] and others [4, 5, 6, 7] that superluminal neutrinos are kinematically 
allowed to emit pairs of e+e~ and would thereby lose most of their energy 
during the CERN-Gran Sasso flight of 730 km [3]. ICARUS experiment [8] 
searched for the Cerenkov emission of electrons in the same CERN-CNGS 
to Gran Sasso neutrino beams, where according to the GC calculation, 63% 
of the neutrinos are expected to decay and no anomalous e^e~ events were 
detected. Using this a much stronger bound 6 < 2.5 x 10~^ can be put on 
the superluminal neutrinos. Another problem with superluminal neutrinos 
which has been discussed [9, 4, 10, 11] is that the change in neutrino energy- 
momentum relation restricts the phase-space for the tt — t- z//i process and the 
pion lifetime would be larger for the OPERA neutrinos produced by the pion 
decay from the CERN-CNGS beam. 

In a measurement by MINOS [12] it was found that muon neutrinos of 
average energy 3 GeV traversing a distance 730 km exceed c by an amount, 
6{E = 3GeV) = (5.1±2.9) x 10^^. This however is in contrast to the neutrino 
observations from supernova SN 1987a [13, 14] where over a flight path of 51 
kpc, the neutrinos with energy in the band (7.5 — 39) MeV all arrived within 
a time span of 12.4 sec and the optical signal arrived after 4 hours of the 
neutrino signal (consistent with prediction of supernova models) from which 
it is inferred that 6{E = IbMeV) < 10~^. 

In the present paper we calculate the rates for e+e~ radiation from neutri- 
nos and pion decay assuming that the neutrinos obey the energy momentum 
relation of the form E"^ = + + ^nP"'- This is motivated by Horava- 
Lifshitz type field theories [15, 16] where the higher derivative terms break 
Lorentz invar iance at high scale but help in removing ultraviolet divergence. 
We find that for n > 2 the bounds on ^„ from direct time of flight measure- 
ment [2] are much more stringent than from constraining the Cerenkov decay 
process during the Cern to Gran-Sasso flight of neutrinos [8]. This is unlike 
the case for n = 2 studied by Glashow-Cohen [3] and [9, 4, 10, 11] where the 
constraint on Lorentz violation from the kinematically forbidden processes is 
more stringent than from the time of flight measurement. 

Models which explain the Opera result of superluminal neutrinos and 
which have a bearing on the question of Cerenkov emission from neutrinos 
or the pion decay kinematics fall broadly in the following categories: 
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1. Deformed Lorentz symmetry models [17, 18, 19, 20, 21, 22, 23] where 
the dispersion relations change from the usual = m? form to a dif- 
ferent form which is still covariant under the modified Lorentz trans- 
formations. In this picture the processes which are forbidden in one 
reference frame (like the rest frame of the massive neutrinos) will be 
forbidden also in the lab frame. 

2. Lorentz invariance violation as in Lifshitz type field theories [15, 16, 24, 
25, 26], from a gauge singlet SUSY sector [27] or a hidden sector [28], 
environmental couplings [29, 30, 31, 32, 33, 34, 35], dynamical symme- 
try breaking [37, 38] , Fermi-point splitting [39] , space-time fluctuations 
[40] and string theory [41]. 

In this paper we we consider the following general dispersion relation 
motivated by Horava-Lifshitz theories 

= + + inP"" (1) 

where n = 2,3,4.. etc. The difference between the superluminal neutrino 
velocity and the speed of light (taken to be 1) is then given by 

5=^-l^^enP"-^ n = 2,3,4... (2) 

The ICARUS time of flight experiment [2] has observed neutrinos and the 
time difference between the neutrino time of flight (tof) and the calculated 
photon tof is 6t = 0.3 ± 4.9{stat) ± 9.3{stat) for neutrino energy E^, = 12.5 
GeV and a distance (731278.0 ± 0.2)m from CERN-CNGS to the detector in 
Gran Sasso. This corresponds to a neutrino superluminality by the amount 

[1], 

6{E = 12.5GeV) = ^^^^ < 2.3 x 10^^ (3) 

From (2) and (3) we put constraints on ^„ from the tof experiment [2]. 

The Cerenkov decay constraint comes from the earlier ICARUS experi- 
ment [8] where the expected number of CC neutrino events is 315 ±5 and the 
observed number is 308. This corresponds to a bound on the decay length 
CT of neutrinos given by 

0.04 < exp(-cr/731.2km). (4) 

We calculate the neutrino Cerenkov decay length cr using the generalised 
dispersion relations (1) and constraint ^„ from the experimental bound (4). 

The upper bounds on the Lorentz violating parameter ^„ for different n 
obtained from the tof experiment [2] and Cerenkov decay experiment [8] are 
displayed in Table 1. 
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n 


Upper bound on ^„ 


Time of flight [2] 


Cerenkov process [8] 


2 


4.6 X 10^^ 


1.7 X 10"^ 


3 


5.9 X GeV-i 


1.3 X 10-6 GeV-i 


4 


3.8 X 10-1° GeV-2 


1.0 X 10"^ GeV-2 


5 


2.5 X 10-11 GeV-3 


7.9 X 10-9 GeV-3 


6 


1.7 X 10-12 GeV-^ 


6.3 X 10-1° GeV"^ 


7 


1.2 X 10-13 GeV"^ 


5.0 X 10-11 GeV-5 



Table 1: Upper bounds on the Lorentz violating parameter ^„ for different n 
obtained from the tof experiment [2] and neutrino Cerenkov decay experiment 
[8] 



2 z/^ ^ z^/iC+e 

We compute the process I'^ij)) — )■ u^{j>')e^{k)e~ {k') for GeV energy neutrinos. 
We will use the dispersion relations (1) for neutrinos in the Lab frame, and 
we will assume that all other particles have the standard energy-momentum 
relations, and we will assume energy momentum conservation in all reference 
frames. This generalises the Glashow-Cohen calculation for the same process 
for the n = 2 energy independent 5 case. 

The amplitude squared for the process is given by 



|M|2 =32^1 [{p-k'){p ■ k) (1 -4sin2^vK + 8sin^^vK)] • (5) 
The decay rate of the neutrino is in general given by 



1 f dy d^k' |M|2 , 

-k) ) (6) 



8(27r)5 7 El E'^ E, 
Without loss of generality we can choose 

p=(E,,0,0,|p|) 
P = {El, Ip'I sin6',0, |p'| cos6') 
k' = El{l, cos sin 6^1, sin0sin6'i, cos6'i). 

The argument of the S function in eq.(6) can be as 

ip~p'- k'f = e„ (p"-^ - (pT"') {p - p) - ppo"" - KD 
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where 

D = in (p""' - - P' + dl) + POI + '^p'OOi COS <P (7) 

Here we will assume that the transverse energy is very small and of the order 
of inP^~'^i therefore we keep only the leading order terms in 9, 9i and 
From here on we shall use the notation p and p' to denote |p| and |p'|. 

Now to fix the limits of the 9"^ and 9\ integrals we need to find their 
maximum values from the 5-function condition i.e 

DE'^ = e„ {p-~' - (pT-') {P - P') - PP'e' (8) 

For the maximum value of 9 we set -Eg = in the above equation so that 
have 

.2 ^ in{p''-'-{pT-'){p-p') .q^ 
"max 

And similarly setting p' = and the electron energy at its maximum i.e 
E'^ = p/2 in the ^-function condition we have, 

{elUa. = ap"-' (10) 

We make the following change of variables to pull out the factors of in 
and p from the integrand : 

p'^Xp, ,9^inP''-^9 ,9^^ir.p''--'9^ 

Using the above definitions in eq. (5) and substituting in eq. (6) gives the rate 
of electron-positron pair emission as 



(^nP""^)^P^ (1 -4sin2^H. + 8sin^^vi/) I dx I ^ d9 

Jo Jo 



d9, / (11) 
Jo 

where / is complicated function of x, ^ and 9i [42]. 

After numerically solving this integral we get the following expression for 
the rate of electron-positron pair emission the general formula for the decay 
width of neutrino splitting process u^p) — > h'{p')e~^{k)e~ {k') comes out to be 

r = ^ (1 -4sin2^H^ + 8sin^^H') J„ [inP''-'f p' (12) 

In eq.(12) /„ is an integral of a function depending on n. The values of /„ 
and r for different values of n are given in Table 2. Using cr for different n 
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n 


2 


3 


4 


5 


6 


7 


8 


In 


1/40 


1/31 


1/29 


1/28 


1/28 


1/28 


1/28 



Table 2: Values of integral in eqn.(12) for different values of n and for 
= 12.5GeV. 

and comparing with the experimental bound from ICARUS [8] shown in (4) 
we obtain the bounds on the Lorentz violating parameter for different n 
displayed in Table 1. 

This generalises our earlier calculation of the n = 2 and n = 4 cases [42] . 
Exact analytical calculations for the n = 2 case has been done [3, 42, 43, 44, 
45]. Our result for the decay width is smaller than the corresponding result 
of [3] and [44] by a factor of 2/3 but is in closer agreement with the results 
of [43, 45]. 



3 Pion Decay 



We calculate the pion decay width in the lab frame with a superluminal 
neutrino in the final state. We assume the dispersion relation E"^ = {p^+inV^) 
in the lab frame. The amplitude squared for the process vr~(g) — )■ fi^ {p)u^{k) 
is, 



IMP = 4G^/X 



ml + ink' 



m 



m 



- + 2 

2 ~ 



(13) 



The decay width is then given by 



8ttE^ 

2 _2 



kdkd cos 6 
q — fcp + 



6{E, + J\q-k\^ + ml-E^) 



< + ink'' 



m 



? + 2 



mf 



(14) 



Writing \q — fcp = k'^ + q^ — 2kqcos6, where 6 is the angle between k 
and q, and E^ = k + ^„A;"~^/2 we see from the argument of the (5-function in 
eq.(14) 

cos d = [ml -ml + 2E^k + ^nfe""'^. - ^fc") {2kq)-^ (15) 
while the derivative of the argument of (5-function with respect to cos 9 yields 



d 



dcos9 



{E, + J\q-k\^ + ml-E, 



kq 



y/k"^ + — 2kq cos 6 + rn?^ 



(16) 
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Pn 


r/Fo 


n = 2 


n = 3 


n = 4 


n = 5 


n = Q 


10 GeV 


0.95 


1.0 


1.0 


1.0 


1.0 


50 GeV 


0.16 


0.48 


0.62 


0.69 


0.73 


100 GeV 


0.04 


0.18 


0.27 


0.32 


0.35 


200 GeV 


0.01 


0.06 


0.11 


0.13 


0.15 


500 GeV 


0.002 


0.02 


0.03 


0.04 


0.05 



Table 3: Ratio of pion decay width in Lorentz violating framework to its SM 
prediction (Fq) at different values of pion momentum (p^) for various n. 



Substituting this in eq.(14) we get 

i + ink- ( 



Glflml f dk 



ml-ml+enk-i'4 + 2 
mr, 



ShEt^ J q 

The limits of the k integral are fixed by taking cos^^ = ±1 in eq.(15) 



;i7) 



kr, 



' 2{E^ - q) 



^"^^^ " 2(E. + q) ^^^^ 

we solve these polynomial equations for k^ax and k^in numerically to obtain 
the kinematically allowed limits of neutrino momentum. Using these limits 
to integrate over the neutrino momentum k we get the decay rate for pion. 
The ratio of pion decay rate thus calculated to the standard model prediction 

ro(,-.,.) = ^(i-^) (19) 

for different n is shown in the Table 3. The n = 4 case has also been dealt 
with in [46] and we are in broad agreement with their result. In Fig. 1 an 
approximate numerical calculation of the pion decay width is plotted as a 
function of pion momentum for different n. 

The decay width for the 100 GeV pions decreases by 65% for n = 6 and 
is smaller at higher n. 
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Figure 1: Ratio of pion decay width in Lorentz violating framework (approx- 
imate numerical calculation) to its Standard Model prediction as a function 
of pion momentum (p^) for different n. 

4 Conclusions 

We have computed neutrino processes assuming a power law correction to the 
neutrino energy-momentum relations. We conclude that for steep power law 
( n > 2) dispersion relations the contraint from the time of flight experiments 
is more stringent than from the measuremnet of the Cerenkov e~^e~ emission 
or from the change in the width of pion decay. Our calculation of neutrino 
Cerenkov emission and pion decay width in Lorentz violating theories can be 
applied for putting bounds on Lorentz violating parameters from the analysis 
of high energy cosmic rays. Also, future experiments for measuring neutrino 
velocities performed at higher energies will put strong constraints on the 
higher derivative Lorentz violation theories [47]. 
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